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Abstract. In the paper we proposed a novel model of unconventional 
computing where a structural part of computation is presented by dy- 
namics of Plasmodium of Physarum polycephalum, a large single cell. 
We sketch a new logical approach combining conventional logic with pro- 
cess calculus to demonstrate how to employ formal methods in design of 
unconventional computing media presented by Physarum polycephalum. 
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1 Introduction 

In the paper we are demonstrating how to design unconventional computing 
media taking into account the problem that in unconventional computing 
media both structural parts and computing data of computers are variable 
ones, in particular in reaction-diffusion processors we are dealing 

with, both the data and the results of the computation are encoded as con- 
centration profiles of the reagents (on the contrary, in conventional models 
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of computation one makes a distinction between the structural part of a 
computer, which is fixed, and the data which are variable and on which the 
computer operates). As a result, while in conventional models circuits are 
fixed, in unconventional computing media circuits could be set up just as 
dynamically variable ones within spatio-temporal logic. Solving this task 
allows us to build up nature-inspired computer models and to consider 
biological and physical systems as computational models. 
One of the unconventional, nature inspired models similar to reaction- 
diffusion computing is chemical machine in that molecules are viewed as 
computational processes supplemented with a minimal reaction kinetics. 
Berry and Boudol first built up a chemical abstract machine f^,'4] as an 
example of how a chemical paradigm of the interactions between molecules 
can be utilized in concurrent computations (in algebraic process calculi). 
We are considering another abstract machine of reaction-diffusion comput- 
ing exemplified by dynamics of Plasmodium of Physarum polycephalum. 
This machine is constructed by using process algebra, too. Using it, we 
are studying a possibility of logical representation of the computation in 
reaction-diffusion systems. Using notions of space-time trajectories of local 
domains of a reaction-diffusion medium we will could the spatial logic of 
trajectories, where well- formed formulas and their truth- values are defined 
in the unconventional way 

Experimental studies and designs of reaction-diffusion computers could be 
traced back to the pioneer discovery of L. Kuhnert (1986). He demon- 
strated that some very basic image transformations can be implemented in 
the light-sensitive Belousov-Zhabotinsky system. The ideas by L. Kuhnert., 
K. L. Agladze, I. Krinsky (1989) on image and planar shape transforma- 
tions in two-dimensional excitable chemical media were further developed 
and modified by N. G. Rambidi (1998). At that time, during the mid and 
late nineties, a range of chemical logical gates were experimentally built in 
the Showalter and Yoshikawa laboratories. The first chemical laboratory 
prototypes of precipitating chemical processors for computational geome- 
try were developed by A. Adamatzky (1996). He also designed and studied 
a range of hexagonal cellular-automaton models of reaction-diffusion ex- 
citable chemical systems. Now reaction-diffusion computing is an extremely 
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Figure 1: An example of computational process in Physarum machine. 
Photographs are taken with time lapse circa 24 hours. 



wide area of researches towards unconventional computing. 

The dynamics of plasmodium of Physarum polycephalum could be regarded 
as one of the natural reaction-diffusion computers. The point is that when 
the Plasmodium is cultivated on a nutrient-rich substrate (agar gel contain- 
ing crushed oat flakes) it exhibits uniform circular growth similar to the 
excitation waves in the excitable Belousov-Zhabotinsky medium (Fig. [T]). 
If the growth substrate lacks nutrients, e.g. the plasmodium is cultivated 
on a non-nutrient and repellent containing gel, a wet filter paper or even 
glass surface localizations emerge and branching patterns become clearly 
visible (Fig. [2| [s]). 

The Plasmodium continues its spreading, reconfiguration and development 
till there are enough nutrients. When the supply of nutrients is over, the 
Plasmodium either switches to fructification state (if level of illumination 
is high enough), when sporangia are produced, or forms sclerotium (encap- 
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Figure 2: Basic components of Physarum and its environment: (a) oat flake, 
(b) propagating pseudopodium, Plasmodium's wave- fragment, (c) oat flake 
colonized by Plasmodium, (d) protoplasmic tube. 



sulates itself in hard membrane), if in darkness. 

The pseudopodium propagates in a manner analogous to the formation of 
wave-fragments in sub-excitable Belousov-Zhabotinsky systems. Starting 
in the initial conditions the Plasmodium exhibits foraging behavior, search- 
ing for sources of nutrients (Fig. [T|. When such sources are located and 
taken over, the Plasmodium forms characteristic veins of protoplasm, which 
contracts periodically. Belousov-Zhabotinsky reaction and plasmodium are 
light-sensitive, which gives us the means to program them. Physarum ex- 
hibits articulated negative phototaxis, Belousov-Zhabotinsky reaction is in- 
hibited by light. Therefore by using masks of illumination one can control 
dynamics of localizations in these media: change a signal's trajectory or 
even stop a signal's propagation, amplify the signal, generate trains of sig- 
nals. Light-sensitive of Plasmodium has been already explored in design of 
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Figure 3: Snapshot of experimental disfi with propagating Plasmodium, 
new activated propagating zone (a) and sites of branching pseudopodia, 
junctions of protoplasmic tubes (b) are shown by arrows. 



robotics controhers g§[2], (|||[3]. 

Despite numerous experimental implementations of Physarum computers 
there is a lack of formalization of the Plasmodium's behavior abstract 
enough to infer high-level principle of information transmission by the Plas- 
modium and accurate enough to reflect peculiarities of the Physarum for- 
aging behavior. In the paper we are trying to fill the gap and offer inter- 
pretation of Physarum behavior in a framework of process calculi. In the 
paper we are trying to define Physarum machine as a process calculus built 
on an unconventional interpretation of logical connectives. 



2 A logical approach to analyzing Physarum machine 

Let us set up the problem how to define logical connectives in Physarum 
machine. The matter is that both structural parts and computing data of 
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Physarum computers are variable. 

Physarum machine may be viewed as a labelled transition system, which 
consists of a collection of states £ = {Pij}i=YN j=rK ^^'^ ^ collection T of 
transitions (processes, actions) over them. Assume T: £ i->- CP(£), where 
= {T: T C £,}. This means that 7{p) consists of all states that a 
reachable from p. The transition system is understood as a triple 



where — >-C £ x T x £ is a transition relation that models how a state 
p E £/ can evolve into another state p' £ L due to an interaction cr G T. 
Usually, {p,a,p') € — > is denoted by p p'. So, a state p' is reachable 



The finite word Q;ia2 • . . a„ is a finite trace of transition system whenever 
there is a finite execution fragment of transition system 

Q = poaipia2 ■ ■ . ctnPn such that pi pj+i for all < i < n. 

The word aia2 ■ ■ - an is denoted by trace{Q). The infinite word aia2 ■ ■ ■ 
is an infinite trace whenever there is an infinite execution fragment of of 
transition system 

Q = poaipia2P2Ct?,Pd, ■ ■ ■ such that pi pi+i for all < i. 
The word aia2 ... is denoted by trace{Q) too. 

Definition 2.1. An infinite (resp. finite) trace of state p denoted by g{p) 
is the trace of an infinite (resp. finite) execution fragment starting in p. 

Each trace can be regarded as a graph, where nodes represent states and 
edges transitions. In this way, transition system is viewed as graph trees. 

Conventionally, logical connectives are defined in the algebraic way that is 
broken within transition systems. Therefore in Physarum machine we will 
distinguish two kinds of logical connectives: 

1. logical connectives defined co-algebraically, these ones are closed to 
conventional logical connectives: they are "eternal," because they are 
defined over traces (i.e. for any future states); 



(£,T,^), 
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2. logical connectives defined as transitions over states, these ones differ 
from conventional logical connectives: they are defined over states, 
therefore their values could change for some future states. 

2.1 Logical connectives defined co-algebraically 

First, let us consider logical connectives defined co-algebraically, i.e. by 
coinduction. 

An infinite trace of state q{p) may be presented as a kind of stream. For a 
trace q{p), we call Qp{0) the initial value of q{p). Wc define the derivative 
gp{0) of a trace q(j>), for all n > 0, by Qp(n) = Qp{n + 1). For any n > 0, 
Qp{n) is called the n-th state of q{p). It can also be expressed in terms of 

higher-order trajectory derivatives, defined, for all A; > 0, by Qp^ = q{p); 
^(^fe+i) _ (^g(^)y_ jj^ ^]jjg gagg ^j^g state of a trace g{p) is given by 

Qp{n) = Q^\o). So, the trajectory is understood as an infinite sequence 
of derivatives: g{p) = gp{0) :: Qpi^) ■'■ £'p(2) :: ••• :: Qp{n — 1) :: gp^\ or 
q{p) = {gp{0),Qp{l),gp{2),...). 

A hisimulation on the set of traces is a relation R such that, for all q{p) and 
Q{q), if {Q{p),Qiq)) G R then (i) ()p{0) = gq{0) (this means that they have 
the same initial value) and (ii) {g'p, g'^) G R (this means that they have the 
same differential equation). 

If there exists a hisimulation relation R with {gp, gq) G R then we write 
gp ~ gq and say that gp and gq are bisimilar. In other words, the bisimilarity 
relation ~ is the greatest bisimulation. In addition, the bisimilarity relation 
is an equivalence relation. 

Theorem 2.1 (Coinduction). For all g{p),g{q), if there exists a bisimula- 
tion relation R with {gip), g{q)) £ R, then g{p) = g{q). □ 

This proof principle is called coinduction. It is a systematic way of proving 
the statement using bisimularity: instead of proving only the single iden- 
tity g{p) = g{q), one computes the greatest bisimulation relation R that 
contains the pair {g{p),g{q)). By coinduction, it follows that g{p) = g{q) 
for all pairs {g{p), g{q)) € R. 
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Now consider logical connectives defined by coinduction over traces. Their 
syntax is as follows: 

Variables: p ::= p \ q \ r 

where p, q, r are states of Physarum machine presented as a labelled tran- 
sition system. 

Constants: c ::= T | _L 

where T means the truth (the ideal, universal trace) and _L means the fal- 
sity (the empty, impossible trace). 

Formulas: ip,ip::=p \ c | \ ip'Vip \ (pAtp \ ip D 

These definitions are coinductive. For instance, 

• a variable p is of the form of a trace p = p(0) :: p(l) :: p(2) 
p(n — 1) :: p*-""-*, where p(i) G {p,q,r, . . .} for each i £ uj; 

• a constant c is of the form of a trace c = c(0) :: c(l) :: c(2) 

c(n — 1) :: c^*^), where c(i) G {T,±} for each i G w, a particular case 
is [T] = [T](0) [T](l) :: [T](2) [T]W, where [T]{i) = T for 

each i E u; 

• a formula -xp has the differential equation {^tp)' = ~'(<^') and its initial 
value is (-i(^)(0) = -'^'(O), this formula will be understood as D [-L]; 

• a formula (p\/ ip has the differential equation {ip V ip)' = ip' \/ ip' and 
its initial value is {^p V ^)(0) = <p{0) V -0(0); 

• a formula ip Aip has the differential equation {ip A ip)' = tp' Atp' and 
its initial value is {(p A ^)(0) = <^(0) A ip{0); 

• a formula ip D ip has the differential equation [(p D ip)' = (p' ip' and 
its initial value is {(p D ip){0) = <^(0) D ip{0)- 
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2.2 Logical connectives defined as transitions over states 

In analyzing the Plasmodium we observe processes of inaction, fusion, and 
choice, which could be interpreted as unconventional (spatial) falsity, con- 
junction and disjunction respectively, denoted by Nil, &; and +. These 
operations differ from conventional ones, because they cannot have a deno- 
tational semantics in the standard way. However, they may be described 
as special transitions over states of Physarum machine: 

1. inaction (Nil) means that pseudopodia has just stopped to behave, 

2. fusion (&) means that two pseudopodia come in contact one with 
another and then merge, 

3. choice (+) means a competition between two pseudopodia in their 
bahaviours. 

Let us notice that a Boolean algebra may be extended up to the case of the 
system of logical connectives defined by coinduction [^[5] (see the previous 
subsection). However, if we define three basic logical connectives (falsity, 
conjunction, disjunction) as transitions over states of Physarum machine, 
they will be extremely non-classical and Boolean properties do not hold for 
them in general case. 



3 Physarum process calculus 

Further, let us try to build up a process calculus combining two approaches 
to logical connectives for describing the dynamics of Physarum machine, 
i.e. we are trying to show that indeed this machine could be presented as a 
labelled transition system with some logical relations. 
Assume that the computational domain is partitioned into computational 
cells Cj = 1, . . . , K such that Cj n Cj = 0, z / j and (jjLi Cj = f^. Then 
suppose that in the K cells, there are active species or growing pseu- 
dopodia and the state of species i in cell j is denoted hy pij, i = 1, ... , N, 
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j = 1, . . . ,K. These states are time dependent and they are changed by 
Plasmodium's active zones interacting with each other and affected by at- 
tractants or repellents. Plasmodium's active zones interact concurrently 
and in a parallel manner. Foraging plasmodium can be represented as a set 
of following abstract entities (Fig. [2]). 

1. The set of active zones (growing pseudopodia or actions) Z = {a, b, . . .} 
(Fig. 2a). On a nutrient-rich substrate plasmodium propagates as a 
typical circular, target, wave, while on the nutrient-poor substrates 
localized wave-fragments are formed. Each action a from Z belongs 
to a state pij, i = 1, . . . , N, j = 1, . . . , i^T of a cell i, which is its 
current position, and says about a transition (propagation) of a state 
Pij to another state of the same or another cell. Part of plasmodium 
feeding on a source of nutrients may not propagate, so its transition 
is nil, but this part can always start moving. 

2. The set of attractants {Ai, A2, . . . } are sources of nutrients, on which 
the Plasmodium feeds. It is still subject of discussion how exactly 
Plasmodium feels presence of attracts, indeed diffusion of some kind 
is involved. Based on our previous experiments we can assume that 
if the whole experimental area is about 8 — 10 ~ cm in diameter then 
the Plasmodium can locate and colonize nearby sources of nutrients. 
Each attract A{a) is a function from a to another action h. 

3. The set of repellents {Ri,R2, . . . }. Plasmodium of Physarum avoids 
light. Thus, domains of high illumination are repellents such that 
each repellent R is characterized by its position and intensity of illu- 
mination, or force of repelling. In other words, each repellent R{a) is 
a function from a to another action b. 

4. The set of protoplasmic tubes {Ci,C2, ■ ■ ■}■ Typically plasmodium 
spans sources of nutrients with protoplasmic tubes/veins (Fig. [2|. 
The Plasmodium builds a planar graph, where nodes are sources of 
nutrients, e.g. oat flakes, and edges are protoplasmic tubes. C{a) 
means a diffusion of growing pseudopodia by an action a. 
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Our process calculus contains the following basic operators: 'NiV (inac- 
tion), (prefix), 'I' (cooperation), '\' (hiding), (reaction/fusion), '+' 
(choice), a (constant or restriction to a stable state), A{-) (attraction), 
R{-) (repelling), C(-) (spreading/diffusion). Let A = {a, 6, . . . } be a set of 
names. With every a G A we associate a complementary action a. Define 
L = {a, a : a G A}, where a is considered as activator and as inhibitor 
for a, be the set of labels built on A (under this interpretation, a = a). 
Suppose that an action a communicates with its complement a to produce 
the internal action r. Define L-j- = LU {r}. 

We use the symbols a, 13, etc., to range over labels (actions), with a = a, and 
the symbols P, Q, etc., to range over processes on states pij, i = 1, . . . , N , 
j = 1, . . . , K . The processes are given by the syntax: 

P,Q::=Nil \ a^P \ A{a)-kP \ R{a)-kP \ C{a) \ {P\Q) \ P\Q \ 
PkQ I P + Q I a 

Each label is a process, but not vice versa. An operational semantics for 
this syntax is defined as follows: 



a-k P — > P 
^ (^(«) = /3), 

- {R{a) = p), 



R{a)irP ^ 

(the conclusion states that the process of the form a-kP (rcsp. A{a) ★P or 
R{a)-kP) may engage in a (resp. A{a) or R{a)) and thereafter they behave 
like P; in the presentations of behaviors as trees, a-k P (resp. A{a) ★ P or 
R{a)'kP) is understood as an edge with two nodes: a (resp. A{a) or R{a)) 
and the first action of P), 
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P P' 

Constant: (a ::= P, a G L,-), 

a — > P' 

p^p> Q^Q' 
Choice: 



P + Q^P' P + Q^Q' 



(these both rules state that a system of the form P+Q saves the transitions 
of its subsystems P and Q), 

P^P' Q^Q' 

Cooperation: 



P\Q ^ P'\Q P\Q ^ P\Q' 

(according to these rules, the cooperation — interleaves the transitions of 
its subsystems), 

P^P' Q^Q' 
P\Q ^ P'\Q' 

(i.e. subsystems may synchronize in the internal action r on complementary 
actions a and a), 

P P' 

Hiding: {a^Q,QCL), 

P\Q^P'\Q 

(this rule allows actions not mentioned in Q to be performed by P\Q), 

Fusion: = — 

a * PkP ^ Nil 

(the fusion of complementary processes are to be performed into the inac- 
tion) , 

pJ^P' Q^P' P^P' Q^P' 

PkQ P' ' QkP P' 

(this means that if we obtain the same result P' that is produced by the 
same action a and evaluates from two different processes P and Q, then P' 
may be obtained by that action a started from the fusion P&iQ or QkP), 



Vol. Logical Modelling of Physarum Polycephalum 13 



P P' P " ) P' 

PkQ ^ Nil + C{a) + P' ' QkP ^ Nil + C{a) + P' 

(these rules state that if the result P' is produced by the action a from the 

processes P, then a fusion PSzQ (or Qk,P) is transformed by that same a 

either into the inaction or diffusion or process P'). 

These are inference rules for basic operations. The ternary relation P 

P' means that the initial action P is capable of engaging in action a and 

then behaving like P' . 

The informal meanings of basic operations are as follows: 

1. Nil, this is the empty process which does nothing. In other words, 
Nil represents the component which is not capable of performing any 
activities: a deadlocked component. 

2. a ★ P, a process a G L followed by the process P: P becomes active 
only after the action a has been performed. An activator a G L 
followed by the process P is interpreted as branching pscudopodia 
into two or more pseudopodia, when the site of branching represents 
newly formed process a ★ P. 

In turn, an inhibitor a G L followed by the process P is annihilating 
protoplasmic strands forming a process at their intersection. 

3. A{a)-kP denotes a process that waits for a value a and then continues 
as P. This means that an attractor A modifies propagation vector of 
action a towards P. Attractants are sources of nutrients. When such 
a source is colonized by Plasmodium the nutrients are exhausted and 
attracts ceases to function: A{a) ★ Nil. 

4. R{a)-kP denotes a process that waits for a value a and then continues 
as P. This means that a repellent R modifies propagation vector of 
action a towards P. Process can be cancelled, or annihilated, by a 
repellent: R(a)-kNil. This happens when propagating localized pseu- 
dopodium a enters the domain of repellent, e.g. illuminated domain, 
and a does not have a chance to divert or split. 
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5. C{a), a diffusion of activator a G L is observed in placing sources of 
nutrients nearby the protoplasmic tubes belonging to a or inactive 
zone (a ::= Nil). More precisely, diffusion generates propagating 
processes wfiicfi establisli a protoplasm vein (the case of activator 
a) or annihilate it (when source of nutrients exhausted, the case of 
inhibitor a). 

6. P\Q, this is a parallel composition (commutative and associative) of 
actions: P and Q are performed in parallel. The parallel composition 
may appear in the case, two more food sources are added to either side 
of the array and then the Plasmodium sends two streams outwards to 
engulf the sources. When the food sources have been engulfed, the 
Plasmodium shifts in position by redistributing its component parts 
to cover the area created by the addition of the two new processes P 
and Q that will already behave in parallel. 

Process P can be split, or multiplied, by two sources of attractants 
(j4i^2)(-P) *-Fi|-F2- Pseudopodium P approaches the site where dis- 
tance to Ai is the same as distance to ^2- Then P subdivides itself 
onto two pseudopodia Pi and P2- Each of the pseudopodia travels 
to its unique source of attractants. Also, process P can be split, or 
multiplied, by a repellent: R{P) *Pi|P2- Biophysics of fission with 
illuminated geometrical shapes is discussed in [5]. The fission hap- 
pens when a propagating pseudopodium 'hits' a repellent. The part 
of pseudopodium most affected by the repellent ceases propagating, 
while two distant parts continue their development. Thus, two sepa- 
rate pseudopodia are formed. 

7. P\Q, this restriction operator allows us to force some of P's actions 
not to occur; all of the actions in the set Q C L are prohibited, i.e. 
the component P\Q behaves as P except that any activities of types 
within the set Q are hidden, meaning that their type is not visible 
outside the component upon completion. 

8. P&lQ, this is the fusion of P and Q; PSzQ represents a system which 
may behave as both component P and Q. For instance, Nil behaves 
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Figure 4: Merging of two Plasmodium's wave-fronts. Photos are made with 
interval 9 hours. 



as P&P, where P is an activator and P an appropriate inhibitor 
respectively. The fusion of P and Q is understood as collision of two 
active zones P and Q. When they collide they fuse and annihilate, 
PSzQ-kNil. Depending on the particular circumstances the new active 
zone a (the result of fusing) may become inactive (Nil), transform to 
protoplasmic tubes (C(a)), or remain active and continue propagation 
in a new direction (the case of prefix 

When two pseudopodia come in contact one with another, they do 
usually merge (Fig. |4]). Thus by directing processes with attractors 
we can merge the processes: A{Pi,P2) *Pi&P2 (see details in [5]). 

9. P + Q, this is the choice between P and Q; P + Q represents a sys- 
tem which may behave either as component P or as Q. Thus the 
first activity to complete identifies one of the components which is 
selected as the component that continues to evolve; the other com- 
ponent is discarded. In Physarum calculi, the choice P + Q between 
processes P and Q sometimes is represented by competition between 
pseudopodia tubes C{P) and C{Q), i.e. C{P,Q) = C{P) + C{Q). In 
other words, two processes P and Q can compete with each, during 
this competition one process 'pulls' protoplasm from another process, 
thus making this another process inactive. The competition happens 
via protoplasmic tube. 
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10. a, constants belonging to labels are components whose meaning is 
given by equations such as a ::= P. Here the constant a is given the 
behaviour of the component P. Constants can be used to describe 
infinite behaviours, via mutually recursive defining equations. 

Thus, in this process calculus we have two kinds of logical connectives. 

1. The group of connectives defined by coinduction. They are derivable 
from the hiding. Indeed, let 1 be a universal set of active zones, then 
the following equalities hold: 

-iP ::= 1\P negation, 

P AQ ::= P\{1\Q) conjunction, 
Py Q::= 1\{{1\P)\Q) disjunction, 
PD Q ::= 1\{P\Q) implication. 
These connectives satisfy all properties of Boolean algebra. 

2. The group of connectives defined as transitions. It consists of three 
operations: inaction, fusion and choice. Their basic properties: 



Nil\P ^ Nil, (3.i; 

PkP ^ Nil, (3.2; 

P&P ^ P, (3.3 

P&iNil ^ Nil, (3.4; 

(P + Q)\P' ^ P\P' + Q\P' , (3.5; 

{PkQ)\P' ^ P\P'kQ\P', (3.6; 

PkQ ^ QkP, (3.7; 

Pk{QkR) ^ {PkQ)hR, (3.8; 

P + P ^ P, (3.9; 
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P + Q ^ Q + P, 
P+{Q + R)^iP + Q)+R, 
Pk{Q + R)^ (PkQ) + (P&P), 
P + (QkR) ^ (P + Q)&(P + R) 



P + Nil^ 



P, 



(3.10) 
(3.11) 
(3.12) 
(3.13) 
(3.14) 



where = is a congruence relation defined on the set of processes. 



4 Conclusion 

In the paper we have just shown that the behavior of Plasmodium of 
Physarum polycephalum could be considered as a kind of process calcu- 
lus with several logical connectives defined in non-standard way. Thus, the 
media of Physarum polycephalum can be viewed as one of the natural un- 
conventional (reaction-diffusion) computers. Its weakest point is that the 
speed of computation is so slow: each new state of Physarum dynamics may 
be observed just in hours (see Fig. [T|. 
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